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Abstract 

A description of a dislocation-free elastic glass phase in terms of domain walls 
is developed and used as the basis of a renormalization group analysis of the 
energetics of dislocation loops added to the system. It is found that even 
after optimizing over possible paths of large dislocation loops, their energy is 
still very likely to be positive when the dislocation core energy is large. This 
implies the existence of an equilibrium elastic glass phase in three dimen- 
sional random field X-Y magnets, and a dislocation free, bond-orientationally 
ordered "Bragg glass" phase of vortices in dirty Type II superconductors. 
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It has been believed for a long time that systems with quenched randomness that couple 
to a continuous symmetry order parameter cannot exhibit long range order in less than 
four dimensions. J!]]- Originally it was argued that positional order of Abrikosov vortex 
lattices in superconductors would be destroyed by random pinning. [I] Recently for a class of 
random systems including X-Y magnets in a random magnetic field, the absence of long 
range order has been proven rigorously. |§ Yet an intriguing open question remains: for weak 
randomness in such systems, are they simply disordered at low temperatures or can phases 
exist which exhibit some kind of topological or other type of order that distinguishes them 
from high temperature disordered phases? This issue has resurfaced in the context of vortices 
in high temperature superconductors; various authors have either implicitly assumed, raised 
the question of, or given qualitative arguments for, the existence of an elastic vortex glass 
phase which is locally lattice like and is free of large dislocation loops. 0- [0 Such a phase 
would probably have power law Bragg-like singularities in its structure factor and in addition, 
have true bond orientational long range order, J7J thus providing a counter-example to the 
general conjecture mentioned above. || 

In the simpler context of three dimensional random field X-Y magnets, Gingras and Huse 
H have explicitly conjectured, and given some numerical evidence in support of, the exis- 
tence of a phase transition to a defect-free phase for weak-randomness at low temperatures. 
Yet at this point, no convincing analytical arguments to support or deny the existence of an 
elastic glass phase have been put forth, |7j , although the delicate balance between elastic 
randomness and dislocation energies has been pointed out by Giamarchi and Le Doussal [|7| . 

In this paper, we explicitly study the stability of a putative elastic glass phase in a three 
dimensional random field X-Y model to dislocation loops — to avoid confusion, we will refer 
to the relevant topological defects in all these systems as "dislocations". To do this we must 
first reconsider the behavior of the ground state and excitations of an elastic glass model 
with dislocations excluded by fiat. The framework that will be developed, naturally allows 
one to analyze the energetics of a dislocation loop that is added to the system. If the core 
energy of the dislocation line is sufficiently large, it is found that, even after optimizing over 
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the possible paths of a dislocation loop, large loops cost energy with high probability. This 
situation can be achieved for a weak random field X-Y model and we thus conclude that an 
elastic glass phase should exist in this system in three dimensions. 

By analogy, our results are applied to other systems, especially the elastic vortex glass. 

Our basic starting point will be the elastic glass model with hamiltonian 

H = g X] - Viv)} 2 -hj^cos [ip(x) - 7O)] (1) 

(xy) % 

with j(x) independent quenched random variables on each site uniformly distributed on 
[— 7r,7r] and f(x)e(— oo, oo). This model has been extensively studied by a variety of tech- 
niques, including an approximate real-space renormalization group (RG) flOU , a perturbative 
4 — e RG expansion |7| and an approximate variational replica symmetry breaking (RSB) 
calculation. |7| At all temperatures the behavior is controlled by a zero temperature fixed 
point whose properties yield disorder averaged (denoted by an overbar) mean square phase 
variations 

< [<p(x) - <p(y)] 2 >a2A\n\x-y\ (2) 

at large distances, with A a universal coefficient computable in the 4 — e expansion. J7| It is 
believed that the mean correlation function will decay as < e lLp ^e^ lip ^ > ~ ( x - y )v (although 
the gaussian variational RSB [jj result that rj — A should not be correct). These results 
have essentially been obtained from coarse graining or Fourier space representation of the 
phase variables. Unfortunately, this framework does not appear to be naturally amenable 
to consideration of dislocations, for these intrinsically involve, as we shall see, physics on 
many length scales. 

A complementary and more complete picture of the elastic glass phase can be constructed 
in terms of domain walls which turn out to be the natural objects at long length scales. 
This is most easily seen by studying the limit h — > oo so that <£>(x) = j(x) + 2irn(x) 
with {n(x)} integers. The ground state can then be represented [up to a uniform shift 
in the {^(2;)}] by the oriented surfaces through which n(x) changes by ±1. The typical 
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<p(x) — f(y) ~ ±JAln(x — y) arises from a sum of random ±2n terms from crossing the 
nested set of closed surfaces enclosing x and y; these occur typically on scales 1, B, B 2 ■ ■ ■ 
with B ~ e 1 /" 4 . In contrast, the coarse grained cp averaged over regions of size of order 
h\x — y\ will vary from x to y by only 0(1). 

We will primarily be interested in studying the ground states of the system with different 
boundary conditions — later including boundary conditions induced by dislocations. Thus 
rather than working with the surfaces across which n(x) changes, it is useful, as for Ising 
spin glasses, |Tl|] to consider configurations relative to the ground state {nc(x)} with some 
chosen fixed boundary conditions; any state can then be represented by the set of oriented 
domain walls across which n(x) — nc{x) changes. A crucial question is the typical energy of 
the minimal domain wall excitation, El, that surrounds a chosen volume L 3 . 

The observation that the coarse grained tp has variation of 0(1) strongly suggests from 
the scaling of the elastic energy in Eq. (|T[), that 

e L ~ L e with 6 = 1 (3) 

(generally 9 = d — 2). This can be seen more explicitly by fixing <p(xi = 0,X2,X3) = in 
a system of size L x L x L and letting <p{x\ = L,X2,x%) change from zero, which defines 
the reference state {nc{x)}, to 2tt, with periodic boundary conditions in X2 — > £2 + L and 
x 3 — y x 3 + L. This forces a single domain wall spanning the system and changes the energy 
by Ei. We can now make use of the powerful statistical symmetry of the model hamiltonian 
Eq. ([!]): if cp is replaced by <p = <pd + i> with ipr, single valued module 2tt and VVd = 
[with lattice derivative operators] then the statistical properties of 

*) = w-^EIW (4) 

with constant boundary conditions on ip are identical to those of TC((p) with constant bound- 
ary conditions on (p. Choosing ^ = ^p-, this implies that the mean energy of the forced 
domain wall is El = |(27r) 2 L. The necessary balance of the random part of the energy 
with the mean elastic part (~ L), implies that this is consistent only if the variations of the 
spanning wall energy, El are also of order L; i.e. 8 = 1. 
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If in the L 3 system, the boundary condition at x\ = L is changed to ip = 2ns, then s 
spanning domain walls will be forced. Because these are closer together for larger s, they 
can less easily optimize their positions and their energy will be larger. Following arguments 
for confined directed polymers, [12|] we consider each section of domain wall of scale L/s. 



Roughly, these are transversely confined by the others on scale L/s. Thus we guess that the 
mean energy of a section will be of order L/s yielding, with s 3 sections, a total increase in 
the mean energy El ~ Ls 2 in agreement with the result from the statistical symmetry. The 
sample-to-sample variations in El can be guessed from a central limit sum of s 3 roughly 
independent variations of each section yielding 5El ~ s(L/s)ys^. For s ~ -|, this yields 
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5El ~ L? as should be expected since the resulting system's ground state — and indeed its 
hamiltonian Ti. — is very different everywhere from the s = reference system. 

The main lesson from the above is that sections of optimal domain walls on scales L 
that are confined on the same scale L, typically have energy that is distributed with positive 
mean and variations both of order L. If we try to put several distinct domain walls of scale 
L into a volume of order L 3 , the energy of each of them will have to increase. 

We now consider inserting a single dislocation loop of radius R into a fixed position of 
the system by making ip multivalued with V x Vtp = 2tt on plaquettes through which the 
dislocation loop passes and V x Vip = elsewhere. If the core energy per unit length of 
the dislocation is e , the statistical tilt symmetry implies that the extra groundstate energy, 
D R , due to the dislocation has mean 

S^=27ri2[e + 7rlni2] (5) 

with the short-scale cutoff of InR absorbed into 6q. How can we understand this in terms of 
domain walls? The dislocation loop forces in a single domain wall that spans the loop. The 
sections of the wall of scale L = 1 adjacent to the loop are confined on this scale and thus 
each have mean energy of order 1 and roughly independent variations of this same order. 
The sections of scale L = 2 are attached to these which confines them on scale L = 2 and 
gives each a mean energy of order 2; and so on, on scales 4, 8, 16 up to L ~ R. Thus each 



factor of two in scale will contribute a factor of R to Dr yielding the RlnR of Eq. (|5|). But 
the variations 5Dr of Dr will be much smaller: from each scale L there will be a random 
contribution ±(i?/L) 5 L, from the sum of order RjL roughly independent variations of 
sections of scale L. The typical variations in Dr are thus 

8D R ~ R (6) 

which are dominated by the largest scale section of the wall. 

For a large fixed dislocation loop, the energy is thus very likely to be large and positive. 
But we must consider the optimization of the dislocation energy over all possible paths 
of the dislocation loop in a region of volume of order R 3 . We do this by an approximate 
renormalization group analysis of the effects of sections of the wall at each length scale on 
the optimal path with minimum energy of a segment of the dislocation loop. 

We focus on transverse deformations of a dislocation segment by distances of order W 
with the effects of smaller scale deformations of the dislocation and the concomitant changes 
in the sections of wall attached to it on scales smaller than W, included in an effective 
mean dislocation energy per unit length iw with local variations around this value. If the 
dislocation were straight on smaller scales, ew would be simply e + vr In W but we expect 
it to be reduced from this by the smaller scale deformations. Our task is to iteratively 
understand how the deformations on a scale W change e on larger scales. 

In the continuum approximation, (valid on large scales) the typical excess energy cost of a 
transverse distortion of a segment of length A of the dislocation by an amount W « A will 
be ~ewW 2 /A from the extra length of the dislocation. Such a distorted dislocation segment 
will have a different spanning wall attached to it than the straight (on scale W) segment. 
But if this spanning wall were completely different over the whole length A, it would imply 
the existence of many (~ A/W) distinct minimal walls in a volume A 3 which is highly 
improbable by the earlier discussion of confined walls. What should be expected, instead, is 
that the minimal spanning walls attached to the two different paths of the dislocation will 
typically only differ in a strip of width W near the dislocation, with all of them attaching 
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to the same optimal wall further away. 

Each section of scale W will be roughly independent and the mean energy of each of 

these sections of the two wall configurations will be the same with differences between them 

also of order W. Thus the total difference in the attached spanning wall energies between the 

i 

energies of the distorted and undistorted segment will be of order ± (J^J 2 W . Balancing 
this difference against the excess effective core energy cost of the distortion from above, 
yields the typical length Aw over which transverse deformations of size W will occur: 

Aw ~ Wei (7) 

which is >> W if ew is large. 

At the next length scale, 2W, the effective e will change. The inclusion of the wall 
segments of scale W increases the mean energy of a length Aw segment by ~ W(Aw/W) 
but the optimization over the scale W deformation of the dislocation decreases it by ~ 
(A w /W)* W. Thus we find that 

i 

hw ~ ew + 7rln2 - a/e w . (8) 

with a some coefficient. 

This is our key result: although the arguments leading to it are not precise, the form 
should be correct for large ew The RG flow of Eq. (|j) implies that the delicate balance |7j 
between the R In R terms in Dr and its ±i? variations, as well as that implied by the almost 
linear growth in Eq. (|7]) of dislocation distortions W\ with A, is resolved by the dominance 
of the deterministic terms in the dislocation energy over even the optimal random ones. 

For sufficiently large eo, the renormalized energy of the lowest energy dislocation loop of 
radius R in a volume ~ R 3 can now be obtained by renormalizing until a scale Wr at which 
Aw ~ R- On longer scales, the optimal dislocation loop of radius R will look essentially 
circular but can still reduce its energy by rotating or moving within the region of volume 
~ R 3 . From scales W in the range R > W > Wr, the renormalization of e will have a 
similar form to Eq. (^) but with the last term replaced by a'(W/R)*. The mean energy of 
the optimal dislocation loop of radius R in a volume of order R 3 will thus be 
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D R ^2nR TrlniZ + eo-O 



e«?,(lni?)^ 



(9) 



for large e . 

The largest scale W ~ i? should dominate the variations in Z}#, yielding 

5D R ~ R«D^. 



(10) 



Large dislocation loops with negative energy will thus be very improbable. [13| Some small 
loops will of course appear for any eo. The effect of these in the presence of an applied 
Vif will decrease the effective long wavelength elastic constant of the vortex glass phase 
and concommitantly the mean energy of large dislocation loops by allowing spanning walls 
with small holes in them. But if eo is greater than some critical value, eo c , these effects 
will only yield finite renormalizations and the elastic glass phase will be stable. In contrast, 
for e < eoc, ew will decrease with length scale eventually becoming negative and leading 
to the proliferation of dislocations of size greater than a correlation length £, even at zero 
temperature. 

So far, we have focused on an unphysical limit of infinite h with eo tuned by hand. 
But the results will apply much more generally. For weak randomness (h « 1) in a three 
dimensional X-Y magnet, the crossover length scale £ p above which the random fields become 
important diverges as £ p ~ h~ 2 . On scales L larger than £ p , the randomness typically 
prevents the system from being able to find more than one optimal configuration in a region: 
another configuration that differed over most of a region of size L 3 would be expected to 
differ in energy by of order h A L 3 which, since this is much larger than the basic energy scale L 
for excitations and stiffness, is highly unlikely. Thus we expect configurations with different 
boundary conditions (or those differing by large dislocation loops) to differ substantially by 
other than multiples of 2ir only on lower dimensional subsets: these will be domain walls 
with thickness ~ £ p . 

Our analysis can be carried through for small h if the system is first coarse grained to 
a scale £ p . This yields an effective core energy e"g p ~ 7iTn£ p from distances within £ p of the 
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dislocation since the energy on these scales is mostly the deterministic elastic energy Larger 
scale deformations of a dislocation are similar to those studied above. If h is sufficiently small, 
€£ p will be large and the energy gained from optimizing the dislocation path will be small; 
the elastic glass phase will thus be stable for weak random fields. 

The static properties of the elastic glass phase are somewhat subtle. The expected power 
law decay of e 1 ^ correlations in the model system should also occur in the elastic glass 
phase of the X-Y random field system. The phase should also exhibit a non-zero elastic 
stiffness: a mean excess energy cost of order L if opposite boundary conditions on cos ip 
are applied at two ends. But its primary distinguishing features will be dynamic as are 
manifested in the analogous elastic vortex glass phase to which we now turn. 

Aside from the complications of dislocations (and domain walls) with different Burgers' 
vectors and of coupling of the real and order parameter spaces, the arguments given above 
for X-Y random field magnets should apply also for three dimensional lattices in weak 
random potentials both for solids in porous media and for vortex lattices in dirty Type II 
superconductors. Power law singularities at Bragg peak positions should be observable in 
the resulting elastic glass phases, especially if the exponent rj is small. Furthermore, the 
absence of large scale dislocation loops means that the bond-orient at ional stiffness will not 
be local and, even though the randomness couples to the local lattice orientation, long range 
hexatic order should exist. [4, 7] But once dislocations proliferate, as the randomness is 
increased, the orientational stiffness will become short range and the randomness coupling 
to local bond orientations will make the system fully disordered. ||] 

The dynamic properties of an ideal dislocation free elastic vortex glass phase have been 
extensively studied. || With only small scale dislocations present, the superconducting 
properties of this phase should persist provided that the phase is "commensurate" with one 



vortex per unit cell; |LlJ if not, the dynamics of the vacancy and the interstitial lines (14 
will dominate the dynamics with the remaining vortices essentially frozen. One might worry 
that dislocation excitations rather than motion of sections of the lattice could dominate 
the dynamics, but their higher energy cost (L In L vs. L) and less likelihood of having low 
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energy, argue against this. 

For weak pinning, a first order transition directly to the commensurate elastic vortex 
glass should occur; this is the residual of the pure system melting transition that will occur 
regardless of the nature of the low temperature phase and is consistent with experiments. 



\ TE \ For stronger pinning, various other phases might exist: a superconducting inelastic 
vortex glass phase with proliferated dislocations; |§ and incommensurate elastic glass phases 
that could be either superconducting or not. Investigation of possible phase diagrams and 
transitions, as well as development of the theoretical understanding of the simpler X-Y 
systems may have to await further progress in methods for analyzing these subtle systems. 

I would like to thank David Huse, Terry Hwa and Matthew Fisher for useful conversa- 
tions. This work is supported in part by the NSF via DMR 9106237, 9630064 and Harvard 
University's MRSEC. 
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